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ABSTRACT 

High-resolution Chandra observations of Tycho's SNR have revealed several sets of quasi-steady, high- 
emissivity, nearly-parallel X-ray stripes in some localized regions of the SNR. These stripes are most likely 
the result of cosmic-ray (CR) generated magnetic turbulence at the SNR blast wave. However, for the amaz- 
ingly regular pattern of these stripes to appear requires the simultaneous action of a number of shock-plasma 
phenomena and is not predicted by most models of magnetic field amplification. A consistent explanation of 
these stripes yields information on the complex nonlinear plasma processes connecting efficient CR accelera- 
tion and magnetic field fluctuations in strong collisionless shocks. The nonlinear diffusive shock acceleration 
(NL-DSA) model described here, which includes magnetic field amplification from a cosmic -ray current driven 
instability, does predict stripes consistent with the synchrotron observations of Tycho's SNR. We argue that the 
local ambient mean magnetic field geometry determines the orientation of the stripes and therefore it can be 
reconstructed with the high resolution X-ray imaging. The estimated maximum energy of the CR protons re- 
sponsible for the stripes is ~ 10'^ eV. Furthermore, the model predicts that a specific X-ray polarization pattern, 
with a polarized fraction ~ 50%, accompanies the stripes, which can be tested with future X-ray polarimeter 
missions. 

Subject headings: ISM: supernova remnants — ^X-rays: individual (Tycho's SNR) — shock waves — turbulence 



1. INTRODUCTION 

X-ray synchrotron emission structures have been observed 
with the superb spatial resolution of the Chandra tele- 
scope in many young s upernova rernnants ( SNRs) (see e.g., 
Vink & Lamind 120031: iBamba et alj l200l_ lUchivama et alj 
5nl2 



2007HPatnaude & Fesenl2009HEriksen et alj20I Ih . The mor- 
phology of the extended, nonthermal, thin filaments observed 
at the SNR edges, and their X-ray brightness profiles, strongly 
support the interpretation that > 10 TeV electrons are ac- 
celerated at the forward shock of the expanding supernova 
shell and produce synchrotron radiation in an amplified mag- 
netic field. Recently, very unusual structures consisting of or- 
dered sets of bright , non-thermal stripes were discovered by 
lEriksen et all (1201 Ih in SN 1572 (Tycho's SNR) with a deep 
Chandra exposure. Understanding these structures presents a 
formidable challenge for current models of X-ray synchrotron 
images of young SNRs. 
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While diffusive shock acceleration (DSA) has long been 
favored as the likely acceleration mechanism for producing 
the highly relati vistic electrons needed for X-ray synchrotron 
emission (e.g., [B landford & Eichler 1987; Jone s & EUisoiil 
II 9911: lMalkov & Drurv 2001; Ellison etal. 200j), a num- 
ber of basic questions remain concerning the nature of DSA 
in SNRs. One is the relative efficiency for accelerating 
protons versus electrons and their dependence on the mag- 
netic field obliquity, another is the origin of the strong 
magnetic field fluctuations required for DSA to produce 
CR ions up to the "knee" near 10'^ eV in SNRs. (e.g.. 
Bell" 1978^. i Blandford & Ei"chieJfT987t fe erezhko & Krvmskii 
20071: IVladimirov et al . 2008; Bvkov et al 



1988; Blasi et al. 



201 lai) . These questions are all interrelated in non-linear dif- 



fusive shock acceleration (NL-DSA) and a consistent model- 
ing of the X-ray stripes seen in Tycho may help resolve these 
problems. 

A fast and efficient mechanism for amplifying magnetic 
field fluctuations with scales below the gyroradii of the en- 
ergetic particles_producing the amplification was proposed 
by iBelll (12004 l2005h . Recently, it has been shown that this 
short-wavelength turbulence can ignite another instability that 
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produces turbulence with wavelengths well above the gyro - 
radii of the responsible CR particles (iBykov et al.l 1201 fbh . 
Large-scale turbulence is particularly important for determin- 
ing the maximum energy CRs that a given s hock can pro- 
duce. Particle-in -cell (PI C) simulations (e.g ., Niemie c et alJ 
I2008t lOhira et al. 2009: .Riauelme & Soitkovskv 2010l) have 
reproduced the basic predictions of IBcI 11 (12004 .' Studies of 
the fast Bell instability accounting for the nonlinear MHD 
cascade of the growi ng rnodes i n a parallel shock precur - 
sor were performed by Bell (2004');'ZirakashvilietalJ |2008|); 
iReville et al. (2008); Zw eibel & E verett ( 2010l) .^ This work 
assumed a CR current as a fixed external parameter and 
showed a fairly broad distribution of the magnetic fluctuations 
due to the nonlinear cascading. The backreaction of the ener- 
getic particle on the shock structure and the turbulence gen- 
eration when a sizable fraction of the shock energy goes into 
relativistic protons is not included in these studies. 

To account for th e CR current back reaction on the fluctu- 
ation amplification, IVladimirov et al.l (|2009) included Bell's 
nonresonant, CR current driven instability in a Monte Carlo 
model of strong shocks undergoing efficient particle acceler- 
ation. The Monte Carlo simulation self-consistently models 
four strongly coupled shock properties: the bulk plasma flow, 
the full particle spectrum, the self-generated MHD turbulence 
including cascading, and particle injection. In a strong shock, 
CRs are produced with enough energy density to modify the 
upstream flow speed. These precursor CRs produce a cur- 
rent that results, via the Bell instability, in the growth of mag- 
netic turbulence. The turbulence, in turn, sets the momentum 
and space dependence of the particle diffusion coefficient and, 
subsequently, the injection and acceleration efficiency of the 
C Rs, closing the sy st em. 

IVladimirov et al.l (l2009l) found that the shock struc- 
ture depends critically on the efficiency of the mag- 
netic turbulence cascading.^ If the cascading along the 
mean field is suppressed (e.g. , [Ooldreich & Sridhar 1992 
iBrandenburg & Nordlundir201 ll) and magnetic field amplifi- 
cation (MFA) is strong, then the CR shock precursor becomes 
stratified, and the turbulence spectrum contains several dis- 
crete peaks well separated in wavenumber, k. We show below 
that these relatively narrow peaks, where the wave spectrum 
energy density kW can be orders o f magnitude greater tha n 
between the peaks (see Figure 3 in IVladimirov etal.1 120091) . 
can produce stripes in synchrotron emission consistent with 
the Tycho SNR observations. The peaks in turbulence energy 
density occur because of the strong coupling between the par- 
ticle acceleration and MFA processes. They will not appear 
in a test-particle calcul ation that ignores this coupling. In 
IVladimirov et alj ( |2009|) . three downstream peaks were found 
for the case with E^n^^ ~ 100 TeV. In other examples without 
cascading, two and four peaks were found with correspond- 
ingly low er and higher Z^max- We note that while the simu- 
lations of IVladimirov et al] (|2009|) were for a parallel shock, 
the Bell instability always produces fastest wave growth along 
the mean field direction, even though the CR current, which 

- In a strong, CR-modified shock precursor, where the magnetic field fluc- 
tuations likely exceed the mean field value, the distinction between a parallel 
and perpendicular shock is blurred. We define a perpendicular shock as one 
where the mean field outside the precursor region is transverse to the local 
shock normal. This mean ambient field outside the precursor determines the 
anisotropy and polarization of the most strongly amplified modes within the 
precursor. Modes with wavevectors along the mean field are most strongly 
amphfied. 

^ Only energy trans fer from long to short wavelengths is included in 
IVladimirov eTalltaOOgh . 



comes fr om the CR density gradient, may be in a different 
direction (lBellll2005h . The shock modification and concentra- 
tion of wave energy into peaks we require should be nearly 
independent of the mean field geometry in all strong shocks 
that show efficient CR acceleration and large amplitude tur- 
bulence as long as cascading is supp ressed. Furthermore, the 
details of the peak structure found bv lVladimirov et aTl ( |2009|) 
are less important for producing strips than the fact that a nar- 
row, long-wavelength peak is generated by the maximum en- 
ergy ions accelerated by the shock. 

A satisfactory explanation of the highly regular, nearly- 
parallel. X-ray stripe s observed in Tycho's SNR by Chandra 
jEriksen et al.l 1201 lb requires a number of conditions to be 
satisfied. The conditions are: (i) the mean magnetic field ge- 
ometry at the outer blast wave, where the stripes are promi- 
nent, must be nearly perpendicular, i.e., the field is perpen- 
dicular to the local shock normal; (ii) the unstable growing 
MHD modes must be linearly polarized and maintain coher- 
ence over a fairly long spatial scale, k, that is somewhat be- 
low the scale size of the CR precursor, (iii) the background 
turbulence must have narrow peaks in wavenumber k; (iv) the 
stripes must persist without significant variation long enough 
to be seen in a deep Chandra observation; (v) the relativis- 
tic electrons producing the synchrotron radiation must have 
a steep spectrum to enhance the emissivity contrast between 
regions of low and high magnetic field; (vi) NL-DSA must 
be efficient enough in a quasi-perpendicular configuration to 
produce a shock precursor structure and MFA; and (vii) tur- 
bulence cascading along the mean field should be suppressed 
to prevent spectral broadening. The model we now detail pro- 
vides all of these properties in rather a natural way. 

2. MODEL 

The coherent nature of the X-ray stripes suggests that 
the underlying magnetic turbulence is strongly anisotropic. 
Isotropic turbulence would not produce extended coherent 
structures with thin stripes. Both the Bell short-wavelength 
instability and the lo ng-wavelength instability, considered by 
iBvkov et al.l ( 1201 Ibh . produce anisotropic turbulence with a 
prominent growth-rate maximum along the mean ambient 
magnetic field direction. The mean ambient magnetic field 
is the field averaged over scales larger than the CR shock pre- 
cursor Since the SNR radius in Tycho is considerably larger 
than the shock precursor, the local ambient field direction may 
vary over the SNR surface. 

The SNR geometry we model is shown in Figure [T] where 
the ambient (i.e., upstream) field. Bo, is in the .x:-direction and 
is tangent to the shock surface for the region of the forward 
shock in the yz -plane. The pattern shown has the linear polar- 
ization vector, b, in the y-direction and the wavevectors of the 
growing modes are along Bo. The rectangular box schemati- 
cally shows the magnetic field amplification region where the 
CR-driven instability occurs. The inflowing fluctuations in 
the ISM field with wavevectors along the mean field direc- 
tion are amplified most rapidly by the CR-current instability 
in the box of scale size l^. This scale determines the coher- 
ence length of the amplified field. The projected scale of the 
box is imprinted into the simulated image shown in the top- 
left inset in Figure [T] and in Figured The linear polarization 
b of the waves propagating along the downstream field of the 
transverse shock results in regions parallel to the shock sur- 
face with alternating high and low field strengths, which will 
translate into bright and dim regions of synchrotron emission. 

In both the Bell and long-wavelength instabilities, the po- 
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Fig. 1 . — Geometry of the simulated supernova shell. The central figure shows half of the shell quarter The upstream magnetic field amplification region, 
where the CR-driven instability occurs, is shown as a box with a scale size that is about the coherence length of the amplified field 1^ . The local mean magnetic 
field in the upstream flow is perpendicular to the flow velocity indicated as the thick arrow. The wavevector of the amplified mode is along the mean field, 
and the hnearly polarized MHD-wave (in the y-direction) has amplitude |b|. The simulated spectrum of the magnetic fluctuations is shown in the left inset. 
We show in the inset the turbulence spectrum, Wik)kRma^/Wo where the normalization factor Wq is expressed through the r.m.s. downstream magnetic field as 
Wo = Sniax (fi^) /(ISti")- The dotted area in the central figure shows a cross-section of the spherical shell filled with the accelerated X-ray radiating electrons. The 
local densities of the Stokes parameters were integrated over the line of sight along the z-axis producing the image illustrated in the right inset and in Fig. [2] 



larization of the turbulence depen ds on the orientation of th e 
ambient magnetic field ( see iBellii2005l: iBvkov et alj|201 Ibl) . 
The amplified field fluctuations have circular polarization for 
a parallel shock where the ambient field and the CR density 
gradient are directed along the shock normal. In contrast, lin- 
early polarized waves, propagating along the ambient field, 
are amplified in the case of a perpendicular shock. The mode 
polarization is important to produce stripes as well as the 
shock geometry. If the stripes were produced in the quasi- 
parallel shock, they would be limited in length to the thick- 
ness of the synchrotron -bright downstream region, AR, and 
this thickness is estimated below to be small compared to the 
length of the observed stripes. The linearly polarized waves in 
a perpendicular shock, however, can produce a nearly static, 
striped emission pattern with a length much longer than AR 
along the shock surface. 

However, even with linearly polarized waves in a uni- 
form perpendicular field, any stripes that might occur will be 
washed out if the waves in the long-wavelength peak have a 
significant spread in wave number, i.e., if Sk > ko, where ko 
is the long-wavelength peak position in the two peak spec- 
trum shown in the inset in Figure [T]"* The peaks found by 

Fluctuations with wavelengths longer than produc ed by the Bell insta- 
bility may also be amplified (e.g., IBvkov et alj|2011bl) or produced by in- 



IVladimirov et alj (|2009|) obey 6k < ko. When combined with 
linearly polarized waves in a perpendicular shock, these nar- 
row peaks result in the creation of stripe-like structures in the 
synchrotron image. 

We constructed synchrotron images of a spherical SNR in 
a model that accounts for the fluctuating magnetic field struc- 
tures described above fol lowing the approach developed in 
IBvkov et all ( l2008i |2009|) . The ampUfied magnetic field in 
the downstream is simulated as a superposition of the pre- 
dominantly linearly polarized modes of random phases, with a 
spectral distribution consisting of a few isolated narrow peaks 
in wavenumber space, as was obtained in the NL-DSA model 
with a CR-driven instability. The magnetic fluctuation spec- 
trum, shown in the inset in Fig. [T] has two peaks; the long- 
wavelength peak at ko IttR'^.^^ (where Rm^x is the gyro- 
radius of the maximum energy ion accelerated by the shock) 
has a width Sk « Q.Sk^. The second peak at the shorter wave- 
length (with ki 33 X ko) is wider. The energy density of 



verse cascading IRiquelme & Sp itkovskv 20lj|). The longer wavelengths 
have slower growth rates then those produced by the fast Bell instability. 
These may somewhat broaden the spectral peak. Any longer wavelength tur- 
bulence will increase the maximum CR proton energy, and if this turbulence 
is included in the Monte Carlo simulation, there would be a shift of kg to 
smaller k. A NL-DSA model that accounts for all of the resonant and nom'es- 
onant instabilities is not yet available. 
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Fig. 2. — Supernova remnant synchrotron emission images simulated in the NL-DSA model, accounting for magnetic field amplification from a CR current 
driven instability. The left panel is the synchrotron X-ray intensity at 5 keV. The degree of polarization of the X-ray emission is shown in the right panel, and the 
values correspond to the color bar. The relatively high polarization fraction is mainly due to the peaked structure of the magnetic fluctuation spectrum and the 
steepness of the distribution of synchi'otron emitting electrons. 



magnetic fluctuations in both peak s are comparable, as wa s 
the case in the NL-DSA model by IVladimirov et alJ (l2009l) . 
The scales corresponding to the short-wavelength peak at ki 
cannot be resolved in Tycho's SNR even with arcsecond res- 
olution X-ray telescopes, while the long-wavelength fluctua- 
tions have scales resolvable with these instruments. 

Regarding the timing and spacing of the stripes, the wave- 
length of the turbulence at A:o must be consistent with the ob- 
served spacing, and the period must be long enough so that 
the structure does not change significantly in a deep Chandra 
observation of duration Afobs- 

In the quasi-perpendicular portion of the SNR shock the 
stripes move transverse to the flow velocity. The structure 
moves with a phase velocity on the order of the Alfven ve- 
locity Vfl and the time scale for changes in the structure is 
^ (A:oVa)"', yielding a constraint (koVa)'^ > Afobs- The sep- 
aration of the stripes is determined by the wavelength ^ l/ko. 
We estimate that the apparent movement of the stripes due 
to this effect is < 0.02"/yr for a SNR at ~ 4kpc, expanding 
in ambient gas of number density 0.1 cm""*, and assuming 
the amplified magnetic field downstream from the shock is 
^ 60 /iG, yielding v,, ^ 450 km s~'. Given the extended spa- 
tial structure of the stripes, their proper motion can hardly be 
noticed in 10 years of Chandra observations. 

An important point to address is why the stripes that appear 
in the X-ray synchrotron emission do not ap pear in the ra- 
dio sy nchrotron observations of Tycho's SNR (iReynoso et al.l 
Il997t IVighet al.llMTTl) . This is naturally explained if the 
> lOTeV electrons producing the X-rays have a steeper dis- 



tribution than the ^ GeV electrons producing the radio. The 
synchrotron emissivity of a power-law electron spectrum with 
spectral index a is proportional to ^("+1 j t^e local emis- 
sivity is relatively very high for large local B and large a. With 
a steep electron distribution, the intensity contrast between the 
peak and trough of the magnetic turbulence produces a greater 
luminosity contrast for the > TeV electrons than for the flatter 
^ GeV electrons. The stripes in radio are not distinct enough 
to be visible in the VLA observations. 

The kinetic model we used here to sim ulate the electron dis - 
tribution is similar to that discussed by iBvkov et alJ (l2009t) . 
The spatially inhomogeneous electron distribution function is 
calculated from the kinetic equation for electrons at a SNR 
shock using a piece-wise parametrization of the particle dif- 
fusion coefficient that is consistent with t hat obtained for ^ 
lOTeV protons by IVladimirov et all(l2009l) . The synchrotron- 
Compton losses of electrons in magnetic fields of y/ (B^) are 
accounted for. The narrow spherical shell of thickness AR 
filled with the accelerated electrons is shown schematically in 
Figure [T] as a dotted region. 

Finally, as IVladimirov et alJ (l2009l) showed, energy cas- 
cading spreads out the peaks in turbulence power and elim- 
inate the key feature needed to produce an ordered pattern of 
stripes. Therefore, such stripes can form only if the turbulent 
cascading along the mean magnetic field is quenched. 

3. DISCUSSION 

The width AR of the region with synchrotron emitting 
electrons in the shock downstream is typically rather narrow 
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(dotted region in Figure \T^. The thickness of the layer de- 
pends on the cooling rate in the downstream magnetic field: 
AR ^ 6 X 10'^ • Vsh8 • fi^G ' ^Tev '^i^' whcrc the r.m.s. mag- 
netic field BniG is in niG, I'shg is the shock velocity in units 
of lO'' km s"', and Zsjev is the electron energy in TeV. Af- 
ter accounting for the projection effect, the value of AR <^ R 
provides the limb brightening observed in many synchrotron 
images of young SNRs. 

While the synchrotron X-ray emission is mainly produced 
in a thin downstream layer that cannot always be resolved in 
X-rays, the transverse (relative to the shock normal) scale size 
of the emitting region is determined by the projected scale size 
of the magnetic field amplification region in the shock precur- 
sor The wave vectors of the fast growing modes flowing into 
the shock must be directed, in the shock rest frame, within 
a fairly narrow cone along the mean ambient field direction. 
The transverse scale size of the locally amplified fluctuations 
should be comparable to the precursor scale size, which is as- 
sumed to be much smaller than the shock radius. It is also im- 
portant that the fast amplification of the field implies that the 
field correlation length /c is also on the order of the precursor 
scale size. Therefore, the length of a stripe is determined by 
the projected scale size of the magnetic field amplification re- 
gion. In the model under consideration, the fast amplification 
of magnetic field fluctuations due to the CR current-driven in- 
stability occurs in the CR precursor of the strong shock. The 
scale size of the CR-precursor, for Bohm diffusion, is a frac- 
tion / of c/vsh X RmsLx, and typically / ~ 0.1 -0.3 from the 
Monte Carlo simulations. 

The distance between the stripes is set by the wavelength of 
the turbulence which is on the order of R^ax- Therefore, the 
ratio of the stripe elongation length to the distance between 
the field maxima is about /c/vsh, and that number also esti- 
mates the number of the stripes in a single coherent set shown 
in Figure |2] 

4. CONCLUSIONS 

It has been known for quite some time that NL-DSA can be 
efficient with a strong coupling between the CR spectrum and 
the magnetic field dynamics. H owever, a nonlinear study , us- 
ing Monte Carlo techniques by IVladimirov et al.l (l2009l) . has 
recently produced a surprising result: if the turbulence cas- 
cading is suppressed along the mean magnetic field, the shock 
precursor becomes stratified, and a multiple-peak structure of 
the amplified, fluctuating magnetic field spectrum forms. 



As we have shown here, the multiple-peak power spec- 
trum, along with the anisotropic character of the CR-current 
driven instability, can result in structures in synchrotron im- 
ages of SNRs similar t o the r ecent observations of Tycho's 
SNR by 'Eri ksen et aH (1201 Ih . The unusual nature of the 
nearly -parallel X-ray stripes is highly constraining and a num- 
ber of conditions must be fulfilled for them to be produced. 
Observationally, the stripe-like structures should be promi- 
nent and resolvable for a quasi-perpendicular shock configu- 
ration, but they are likely too small or indistinct to be resolved 
with current instruments in regions where the shock is quasi- 
parallel. The synchrotron emission from the stripes should 
be linearly polarized. The simulated degree of polarization is 
shown in the right panel in Figure |2] The high polarized frac- 
tion of ^ 50% predicted for the bright stripes should make 
these structures highly prominent in future X-ray polarization 
observations. The lack of Faraday rotation and the relatively 
steep spectra of synchrotron emission are adva ntageous for 
X-ray polarization studies of young SNRs (see iBykov et al.l 

llool 



Stripe-like structures should form in a section where the lo- 
cal field lies along the shock surface and where the turbulence 
cascading is suppressed. We argue that it is difficult to pro- 
duce the stripe structure in synchrotron images in any natural 
way other than with narrow peaks in the magnetic turbulence 
in a perpendicular shock. This uniqueness offers a new way 
to infer the geometry of the local ambient magnetic field. The 
magnetic field structure is particularly important in modeling 
the TeV gamma-ray emission from Tycho's SNR recently de- 
tected with the VERITAS g round based gamma-ray observa- 
tory bv .Acciarietal] (|20I lb . 
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